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Abstract 

In this work we use two different but complementary approaches in order to study the ghost 
propagator of a pure SU(3) Yang-Mills theory quantized in the linear covariant gauges, focusing 
on its dependence on the gauge-fixing parameter ^ in the deep infrared. In particular, we first 
solve the Schwinger-Dyson equation that governs the dynamics of the ghost propagator, using a 
set of simplifying approximations, and under the crucial assumption that the gluon propagators 
for ^ > 0 are infrared finite, as is the case in the Landau gauge =0). Then we appeal to the 
Nielsen identities, and express the derivative of the ghost propagator with respect to ^ in terms of 
certain auxiliary Green’s functions, which are subsequently computed under the same assumptions 
as before. Within both formalisms we find that for > 0 the ghost dressing function approaches 
zero in the deep infrared, in sharp contrast to what happens in the Landau gauge, where it known 
to saturate at a finite (non-vanishing) value. The Nielsen identities are then extended to the case 
of the gluon propagator, and the ^-dependence of the corresponding gluon masses is derived using 
as input the results obtained in the previous steps. The result turns out to be logarithmically 
divergent in the deep infrared; the compatibility of this behavior with the basic assumption of a 
finite gluon propagator is discussed, and a specific Ansatz is put forth, which readily reconciles 
both features. 

PACS numbers: 12.38.Aw, 12.38.Lg, 14.70.Dj 
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I. INTRODUCTION 


The infrared (IR) behavior of Yang-Mills Green’s functions in the Landau gauge has 
been the subject of numerous studies in the past few years, both in the continuum and 
on the lattice. Particularly important in this challenging quest has been the two-point 
sector of the theory, where it has been hrmly established [1-5] that the gluon propagator 
saturates in the deep IR, a behavior directly associated with the dynamical generation of a 
momentum-dependent gluon mass [6-13], and that the ghost propagator remains massless, 
being accompanied by a dressing function that reaches a hnite value at the origin [14, 15]^. 
Interestingly enough, these characteristic features persist when implementing the transition 
from pure Yang-Mills to real world QCD; specihcally, the inclusion of a small number of 
dynamical light quarks induces quantitative but not qualitative changes to the gluon and 
ghost propagators [20-24], 

Given that the Green’s functions depend on both the gauge-hxing scheme employed and 
the choice of the gauge fixing parameter (gfp), it is important to explore their main dynam¬ 
ical features in different gauges, in order to hlter out the truly gauge-independent properties 
of the theory. In particular, it would be interesting to establish the extent of validity and the 
possible modifications induced to the underlying mechanisms that endow the fundamental 
degrees of freedom, namely quarks and gluons, with their corresponding dynamical masses. 
Furthermore, even though physical observables are ostensibly gauge-independent, nonper- 
turbative calculations are subject to truncations, which in turn may distort the delicate 
conspiracy of terms that produce the required gauge cancellations. It would be therefore a 
useful exercise to probe explicitly the gauge-(in)dependence of certain special combinations 
of Green’s functions that are extensively used in a variety of phenomenological applica¬ 
tions [25-31]. 

Among the different classes of gauges, the linear covariant (or R^) gauges [32] hold a 
prominent position. The corresponding gauge-hxing term that must be added to the stan¬ 
dard Yang-Mills Lagrangian is given by where ^ represents the gfp; some char¬ 

acteristic values include the aforementioned Landau gauge (.^ = 0) and the Feynman gauge 
(.^ = 1). gauges have the advantage of manifest Lorentz covariance, and are particularly 

^ For additional studies and alternative approaches, see e.g., [16-19] and references therein. 
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easy to use in diagrammatic calculations. In addition, by using the novel algorithm proposed 
in [33], they can be implemented in numerical simulations of lattice regularized Yang-Mills 
theories even for 7 ^ 0 [34], 

In the present work we initiate a study of the IR dynamics of the Yang-Mills two-point 
functions within this latter class of gauges, with the main objective to go beyond the standard 
Landau gauge paradigm. To that end, we will resort to two distinct but complementary 
approaches: on the one hand the Schwinger-Dyson equations (SDEs) [35] of the theory, and 
on the other the so-called Nielsen identities (NIs) [36, 37]. 

Within the SDE context, we focus exclusively on the integral equation governing the 
dynamics of the ghost dressing function, F{q^), which has a much simpler structure than 
the corresponding equation for the gluon propagator.At the formal level, the SDE in question 
is written down for general and after approximating the ghost-gluon vertex by its tree-level 
value, the solutions are obtained for the range 0 < .^ < 1 , thus spanning the values between 
the Landau and the Feynman gauges.Our main Ending is that, contrary to what occurs in 
the Landau gauge, F{q^) vanishes as —)■ 0 for all values of ^ within the aforementioned 
interval. This drastic change in the infrared behaviour of F{q^) away from the Landau gauge 
may be traced back to the massless contributions associated with the .^-dependent part of 
the gluon propagator entering into the ghost SDE. Specifically, even if one assumes that the 
cofactor A(g^) of the transverse part of the gluon propagator is finite in the deep IR (as 
happens in the Landau gauge), it is a text-book fact that the longitudinal part (proportional 
to receives no quantum corrections, and maintains its tree-level form [see Eq. (2.2)]. This 
massless contribution, in turn, introduces an infrared divergence into the ghost SDE, which, 
within the approximations employed, can be counteracted only if the solution for F{q^) 
vanishes in the deep IR. In particular, as we will see in detail, F{q^) vanishes at the very 
mild rate of (—c.^logg^//i^)“^/^ (with c > 0 ). 

We then turn to the NIs, which express the gauge-dependence of ordinary Green’s func¬ 
tions (propagators, vertices, etc.) in terms of special auxiliary functions associated with the 
extended Becchi-Rouet-Stora-Tyutin (BRST) sector of the theory^. In the case of the ghost 

^ The gfp-dependence of Green’s functions can be in principle obtained also by using the so-called Landau- 
Khalatnikov-Fradkin (LKF) transformations [38, 39]. These transformations have been used only in an 
Abelian context and are in general formulated in position space; therefore, their use for the problem at 
hand appears to be less direct. 
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dressing function, the corresponding NI permits us to estimate its first derivative of F{q^) 
with respect to for arbitrary values of however, for practical purposes we limit our 
analysis to those ^ that satisfy the condition 1. The reason for this choice is that, in 
this particular limit, the auxiliary functions appearing in the NI may be computed in their 
one-loop dressed approximation, using as input the gluon and ghost propagators known 
from the Landau gauge. The emerging expressions, when evaluated in the deep infrared, 
reproduce rather faithfully the behavior obtained from the ghost SDE; specihcally, up to a 
multiplicative factor, one recovers precisely the derivative of (—c.^logg^//i^)“^/^ with respect 
to i- 

Finally, taking advantage of the Nl-based machinery developed here, we go one step 
further, and study the .^-dependence of the gluon two-point function, which, in the low 
momentum region under scrutiny translates directly into a statement on the dynamically 
generated gluon mass. The relevant auxiliary functions are evaluated using again the ap¬ 
proximations and assumptions employed in the previous case. The result reveals that the 
.^-derivative of the gluon mass displays an IR logarithmic divergence, which can be traced 
back to the masslessness of the ghost propagator. As we explain in terms of an explicit ex¬ 
ample, such a divergent derivative may originate from perfectly IR hnite gluon propagators, 
such as those found in the lattice simulations of [34] for 1. 

The article is organized as follows. In Sect. II we set up the ghost gap equation, discuss 
the approximations and assumptions employed, and present its numerical solutions, paying 
particular attention to the deep IR behavior. In Sect. Ill we address the same problem from 
the point of view of the NIs. Focusing on the identity satished by the ghost dressing function, 
we evaluate it numerically within the one-loop dressed approximation, which allows for the 
determination of the leading IR behavior of F. The result turns out to be in excellent 
qualitative agreement with that found in the previous section. In Sect. IIIC the NI analysis 
is extended to the gluon propagator. In particular, a constraint on the IR behavior of the 
dynamical gluon mass is obtained, and an Ansatz for the possible .^-dependence of the gluon 
mass is proposed. Our conclusions are presented in Sect. IV. Finally, the technical details 
necessary to derive the Yang-Mills NIs are summarized in Appendix A. 
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II. SCHWINGER-DYSON EQUATION ANALYSIS 

In this section we carry out a general analysis of the SDE that governs the ghost propa¬ 
gator, and eventually its dressing function. 

A. General considerations and approximations 

The ghost gap equation (Fig. 1) can be obtained directly from the one-loop ghost self¬ 
energy equation by fully dressing the internal gluon and ghost lines and one of the gluon ghost 
vertices appearing in it [35]. Dressing the right vertex, the SDE for the ghost propagator in 
a linear covariant gauge reads (factoring out the trivial color structure 5“^) 

D-\q^) = 

= q^ + ig'^CA [{k + qYD{k + g)A^^(/c)r''(/c -h g, -k, -q), (2.1) 

Jk 

where n(g^) represent the ghost self-energy, Ca is the Casimir eigenvalue of the adjoint 
representation, and the integral measure is dehned as ^"^k, with jj, the ’t 

Hooft mass and d = 4 — e the dimension of the space-time. and D denote, respectively, 
the gluon and ghost propagators, dehned according to^ 

iA^^(g) = -i P^^(g)A(g2) + ; Ptiu{q) = 9t,u - 

( 2 . 2 ) 

qz 

where S, is the non-negative gfp [32] (see also Appendix A 1), and P(g^) is the so-called ghost 
“dressing function”. E*^ represents the full ghost-gluon vertex, with (all momenta entering) 

r''(gi + g2, -gi, -92) = A{qi + g2, -gi, -g2)g2 + ^(di + d2, -gi, -g2)gi, ( 2 . 3 ) 

where gi (g 2 ) is the gluon (antighost) momentum; at tree-level, = 1 and = 0. 
Then, using Eqs. (2.2) and (2.3), we may rewrite Eq. (2.1) as 

= f + ig^CAq^q'' I D{k + q)A{k)Pi,^{k)A 

+ <q^-CAlD(k + q)(l + Y)[B+Y-Pj, ( 2 . 4 ) 

^ Our conventions can be found in Appendix A. 
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FIG. 1: The ghost gap equation. White (respectively, black) blobs represent connected (respec¬ 
tively, one-particle irreducible) Green’s functions. 

where the common argument {k+q, —k, —q) of the form factors A and B has been suppressed. 

Solving this equation in its full generality would require either independent knowledge 
of the gluon propagator and the form factors of the ghost vertex for general or to couple 
(2.4) to the corresponding SDEs describing A, A and B. However, apart from the lattice 
study of [34], which investigated the gluon propagator for very small values of .^ (.^ < 10“^) , 
there is no direct knowledge of the aforementioned quantities. As for solving the full coupled 
system of SDEs, unfortunately it constitutes a task that lies beyond our present powers. 

Therefore, we will instead study the SDE of Eq. (2.4) within the one-loop dressed approx¬ 
imation, which is obtained by keeping the propagators fully dressed and assigning tree-level 
values to A and B. In addition, we will approximate the A(g^) appearing in the hrst term 
on the right-hand side (rhs) of Eq. (2.4) by the Landau gauge propagator AL(g^). The main 
underlying assumptions behind this later approximation are that A(g^) saturates in the IR, 
assuming the standard form 

A-\q^) = qM{q^)-m^{q^), (2.5) 

and that the deviation between A“^(g^) and AL(g^) in the intermediate momenta region 
is relatively mild, at least for 0 < .^ < 1. Of course, as one approaches the region of 
larger momenta, the perturbative behavior will eventually set in; at one-loop order, A“^(g^) 
renormalized in the momentum-subtraction (MOM) scheme is given by 

A"^(g2) ~ gV(g^) = 1 + log ^ , (2.6) 

where fi is the renormalization point. For example, for the typical values used in this work, 
i.e., jj, = 4.3 and = 0.22, the difference between the Landau and Feynman gauge 

perturbative tails is no more than 7% in the momenta range 2 -=- 5 GeV. In any case, as will 
become clear in the ensuing analysis, the behavior of F{q^) in the deep IR is not particularly 
sensitive to the above considerations; in fact, the complete knowledge of the gluon propagator 
would only affect the subleading terms. 
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Thus, the simplihed version (2.4) that we will consider is given by 


=q^ + ^ D{k + q) q‘‘q''A,(k)P„,(k) + { tl + ^ 


This particular integral equation must be properly renormalized, through the introduction 
of the appropriate renormalization constants for D, Al, and As is well-known, in prin¬ 
ciple the complete renormalization procedure must be carried out multiplicatively. As a 
result, in addition to the ghost renormalization constant Zc that will multiply the tree level 
term q^, further constants multiplying the remaining terms on the rhs of Eq. (2.7) must 
be included; this, in turn, adds an inordinate amount of complexity to the entire prob¬ 
lem. Following the standard approximation, we will simply replace —)■ Zc,q^, and set all 
multiplicative constants equal to unity, thus employing subtractive instead of multiplicative 
renormalization [40, 41]. The actual expression for Z^ is hxed from Eq. (2.7) through the 
momentum subtraction (MOM) renormalization condition Zl“^(/i|) = /i|, where /i| is the 
renormalization point. 

As an elementary check, we may recover from Eq. (2.7) the one-loop expression for F{q^). 
In particular, setting tree-level values for D{k + q) and Al(A;), it is straightforward to show 
that 

F-\q^) = 1 + [(3 - 0 / , 2 .,^ ^2 + 2(1 - 0 / • ( 2 - 8 ) 

4 ^ Jkk\k + qf Jkk\k + qf\ 

Using standard integration formulas, setting and renormalizing in the aforemen¬ 

tioned scheme, one obtains for the renormalized ghost dressing function 


^R^(gE) = i + ^(3-0iog(g^//i^), 


where we have dehned o;* = g^/47r. 


B. Numerical analysis 


After a set of basic algebraic manipulations, together with the shift k + q ^ k, we may 


cast Eq. (2.7) in the form 


( 2 . 10 ) 


This last form of the ghost SDE is more convenient for the numerical analysis that follows, 
because it allows us to carry out exactly the angular integration in the term proportional 
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to while in the hrst term the angular dependence has been passed from the unknown 
function D{k + q) to the function Al(A;), which is known from the lattice. 

In order to solve this equation, we hrst pass to Euclidean space using the standard 
substitution rules 


-»idX; (-ql, -kl,-k^ ■ fc); A((/^), D(q‘‘) -Ae(«S, -Be(9e). 

( 2 . 11 ) 

and suppress throughout the subscript “E” in what follows. Next, we introduce spherical 
coordinates (in d = 4), through the relations 


X = q ] y = z = {k + q) = a; + ?/+ cosd; 


'/Ce 


( 27 r) 


ddsin^d / dyy, 


( 2 . 12 ) 


use the result 


de 


sin^ 6 71 


-Q{x-y) + 


-Q{y-x) 


(2.13) 


Jo 2 [x{x-y) yiy-x) 

where Q{x) is the Heaviside function, and factor out a from both sides of Eq. (2.10). Thus, 
we obtain the hnal equation for the (subtractively renormalized) ghost dressing function 
F{x), 

(XsCa 


F-\x) = Z,- 

. A 


^ sin^ 0 


+ ?- 


IOtt 


Jo 


dyyF{y) + dy 


£M 

V 


(2.14) 


Before proceeding to the full numerical treatment of this integral equation, it would 
be useful to identify some of its main IR features by means of a more direct method. In 
particular, if we assume that the F{x) reaches a hnite value in the IR (a; —)■ 0), inspection of 
Eq. (2.14) reveals that the dominant term in that momentum region is the last one. Indeed, 
the hrst term corresponds qualitatively to the Landau gauge case: if the gluon propagator 
(Al) saturates in the IR, this term is hnite. The second term is also hnite in the IR, as the 
simple change of variable y = tx immediately demonstrates. Therefore, keeping only the 
dominant IR contribution on the rhs of Eq. (2.14) we obtain 


^ ^ X^O^ 1671 L ^ 


y 


(2.15) 




FIG. 2: (color online). The lattice SU(3) gluon propagator evaluated in the Landau gauge [4] and 
the corresponding fit used in our calculation [27]. The dashed curve shows a fit featuring an IR 
maximum which is due to the presence of (divergent) contributions to the gluon (inverse) dressing 
function [42]. All functions are renormalized at /r = 4.3 GeV. 


This integral equation can be converted into a differential equation, by differentiating 
both sides with respect to x; we then obtain 


F’(x) ~ (c 

x —>-0 


F^{x) 


c = 


X 


OIsCa 

IbTT 


( 2 . 16 ) 


which is solved by 

F{x) ~ ± ^ (2.17) 

^Ja — 2^clog(x//i2) 

with a a (possibly ^ dependent) constant, and /i a suitable renormalization scale; the physical 
solution corresponds to the positive sign. Notice that the IR solution given in Eq. (2.17) 
requires the aforementioned non-negativity condition ^ > 0, since otherwise F would become 
complex; in particular, from now on, we will restrict our attention to G [0,1]. 

Eq. (2.17) predicts an important qualitative modification in the IR behavior of the ghost 
dressing functions, compared to what is known from the Landau gauge studies. Specihcally, 
whereas in the Landau gauge El( 0) =const, whenever > 0 one hnds that F is driven to 
zero at the origin, namely F{0) = 0. 

We next focus on the complete numerical evaluation of Eq. (2.14). To this end, we will 
use as input for Al the fit to the available SU(3) lattice data [4] introduced in [27] (see 
Fig. 2). The value of the renormalization point within the MOM scheme is /i = 4.3 GeV. 
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FIG. 3: (color online). Solntion of the SDE (2.14) (upper panel) and the associated ghost prop¬ 
agator (lower panel) for various values of the gauge fixing parameter In the IR the solution 
obtained is perfectly described by Eq. (2.17) after fitting for determining the value of the arbitrary 
constant a. Eor comparison we plot also the Landau gauge lattice data of [4]. 

Notice that in Fig. 2 we show also a fit displaying an IR maximum that must appear due 
to the presence of divergent terms contributing to the gluon (inverse) dressing function [42] 
(see also Sect. IIIB); however, the results finally obtained from the solution of the SDE 
are completely insensitive to the implementation of this particular feature in the gluon 
propagator. 

The solutions obtained for ag = 0.29 and gfp values ranging from 0 to 1 are shown in 
the left panel of Fig. 3. The value of is chosen so that in the Landau gauge ^ = 0 


10 





















one reproduces the lattice data of [4] (see the black continuous curve in Fig. 3); the 30% 
deviation from the expected value of as = 0.22 (at fi = 4.3 GeV) is due to the use of the 
tree-level ghost-vertex, as demonstrated in [43]. 

One immediately observes the drastic change in the IR behavior of the ghost dressing 
function: at .^ = 0 Rl( 0) is hnite, whereas when % 0 F(0) vanishes. The IR behavior is 
precisely the one described by the IR solution (2.17), where 

a = a(0 = 0.12(1+ 0; c = 0.035. (2.18) 

Evidently, the rate at which F{q^) approaches zero is very slow, and begins to set on at the 
rather low scale of about 100 MeV (upper panel of Fig. 3). However, the first appreciable 
deviations from the obtained in the Landau gauge manifest themselves at the higher 

scale of about 300 MeV, where the F{q^) displays a characteristic maximum. This particular 
feature, in turn, may serve as a guiding signal in future lattice simulations away from the 
Landau gauge. 

The overall effect of F{q^) on the full ghost propagator D{q^) is shown in the right panel 
of Fig. 3. In particular, one observes that the rate of divergence of the ghost propagator at 
the origin becomes slightly softer compared to that of the Landau gauge. 

Let us conclude this section by determining for later convenience the IR behavior of the 
derivative with respect to ^ of the ghost dressing evaluated at = 0; one hnds 

~ CsDslog— X Rl( 0 ); Csde = —TT ( 2 . 19 ) 

^ IbTT a(U) 

where we have used the fact that Rl( 0) = l/\/a(0). Clearly, this quantity displays an 
IR logarithmic divergence; substituting the numerical values of the constants involved one 
obtains Csde = 0.15. 

III. NIELSEN IDENTITIES 

In this section we take a different but complementary look at the problem, by resorting 
to a set of identities originally introduced by Nielsen [36, 37]; for all technical details the 
reader is referred to Appendix A, where the general derivation is summarized. 
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A. Ghost propagator 


Consider the ghost two-point sector of the theory. The corresponding NI is readily ob¬ 
tained by differentiating the functional identity (A20) with respect to one antighost and one 
ghost held; setting afterwards all helds and sources to zero, one obtains the relation 


where 


d^V^a^b{q^) = -q)V^d^b{q^), (3.1) 


= -i5^^q^F-\q^)-, = -Yi{q^). (3.2) 


In Eq. (3.1) 0* denotes the antiheld associated to the held 0. In addition, y represents the 
static (he., momentum independent) source associated to the gfp therefore, and despite 
their appearance, all functions in the identity above are two-point functions. 

Eq. (3.1) can be further simplihed by noticing that the so-called ghost (or Faddeev-Popov) 
equation (A17) yields 

, n 

(3.3) 




a result which allows to trade the function P^a* in (3.1) for a ghost two-point function Peg- 
Then, factoring out the trivial color structure 5“^, one is left with the identity 


5?reg(g2) = 




^rg^A^(g,0,-g) +iPexc*(g,0,-g) 


q 




(3.4) 


In order to appreciate with a concrete example how the NIs work, let us consider the 
explicit realization of Eq. (3.4) at the one-loop level. The left-hand side (Ihs) of Eq. (3.4) 
can be immediately deduced from Eq. (2.8), yielding 

1 _ 7 k-q 




g^c 


A 2 

q 


Uk k^{k + g)2 ^ ^ Jk k\k + g)^ J 


(3.5) 


Turning to the rhs of (3.4), the diagrams contributing to the auxiliary functions Pg^A^ and 
at one-loop level are shown in Fig. 4. Using the Feynman rules reported in Appendix A 


exC 


and Ref. [44], one has the results 






ir^^c*{q,^,-q) = * 


.g^CA f k^ + k-q 


Ik k^{k + qy' 


(3.6) 
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FIG. 4: One-loop diagrams contributing to the auxiliary functions and Fc^c* appearing in 

the ghost two-point Nielsen identity (3.1). Notice the presence of the mixed propagator l^bA- 

Notice that the contribution proportional to ^ that could be in principle gener¬ 
ated from diagram (6) of Fig. 4 vanishes as a result of the Slavnov-Taylor identity 
= 0. Thus one dually has 


which, in view of Eq. (3.5), confirms the validity of Eq. (3.4) at one-loop. 


k-q 


Ik k*{k + qy 

(3)7) 


? r(i) 


q^ 




■p(o)/ 2\ 9 Ca 2 

ry(g ) = —^q 


Ik k^ik + qy 


+ 2 


B. Small ^ limit and the one-loop dressed approximation 


Consider now the limit 1; in this case, one can set .^ = 0 on both sides of Eq. (3.4), 
and use Eq. (3.2) to obtain 


d^F{qy = 


^'^5xaS 9,0, -q) + iT^cxc*i9,0, -q) 


F^q^ 


(3.8) 


where the auxiliary ghost functions appearing on the rhs are now evaluated in the Landau 
gauge (see also the discussion at the end of Appendix A 1). This last equation can be used 
to deduce the IR behavior of E(g'^,0 from the knowledge of the basic Green’s functions in 
the Landau gauge. In particular, it allows us to compare the result obtained from the direct 
evaluation of the rhs of Eq. (3.8) in the limit —)■ 0 with the corresponding expression 
derived in Eq. (2.19) in the SDE context. To this end, we will study the auxiliary functions 
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and F^^^* in the one-loop dressed approximation, in which the diagrams contributing 
to each function are obtained from those shown in Fig. 4 by fully dressing the propagators, 
while keeping all vertices at their tree-level values^. The simple inspection of the diagrams 
given in Fig. 4 suggests that, indeed, a logarithmic behavior similar to that of Eq. (3.8) is 
expected to make its appearance. This is because diagrams (a) and (c) may be essentially 
regarded as closed ghost loops, which, due to the nonperturbative masslessness of the ghost 
propagators entering in them, are known to diverge logarithmically in the IR [42, 45]. 

Let us then evaluate explicitly the one-loop dressed expressions of Fg^^^ and one 

has the following results 


.q'^Ca \ f {k-q){k-q+q‘^) 




■Fi^{k)A-^{k +q) , 


— {k-qy 




The terms proportional to the product of two ghost dressing functions Fl in both functions 
are those corresponding to the aforementioned ghost-loops; therefore, in the deep IR both 
functions display a logarithmic divergence, so that, in turn, one has 

r CNilog^ X Fl(0), (3.10) 


where Cni a suitable constant and /i the renormalization scale chosen. Notice that this is 
exactly the kind of behavior found in Eq. (2.19) from the SDE analysis. 

The qualitative agreement between Eq. (2.19) and Eq. (3.10) motivates a further quan¬ 
titative study, focusing on the actual value of the coefficient c obtained within the two 
methods (SDE vs NI). To accomplish this, we evaluate numerically the one-loop dressed 


^ Note that the 6-equatioii (A14) implies that every Green’s function which involves the Nakanishy-Lautrup 
multiplier b remains fixed at its tree-level value: therefore in the Fsector the one-loop dressed approxima¬ 
tion is exact. 
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FIG. 5; (color online). Contributions of the one-loop dressed auxiliary functions to the ghost two- 
point function Nielsen identity. The IR region is perfectly described by the predicted Cni log 
behavior yielding c^i = 0.33. 

contributions (3.9), which are given by (Euclidean space) 

iS, I " + \/f) i^x(!/)n(.-) 

piT roo 

+ fg / d^sin^^/ dyF^{y)Aj^{z) = {a) + {b), 

3(27r)'^ Jo Jo 

d» (l + ^cosej-^FMFM-{c), 

( 3 . 11 ) 

where (a), (6) and (c) denote the contributions of the diagrams appearing in Fig. 4. At 
this point all integrals can be evaluated provided that we supply as input the Landau gauge 
gluon propagator Al and the ghost dressing function Fl (see Fig. 2 and Fig. 3, respectively). 

The results obtained for the three individual terms (a), (6) and (c) of Eq. (3.11), as 
well as their sum, are shown on the left-panel of Fig. 5. One sees that terms (a) and (c) 
show the claimed logarithmic divergence, while in the case of (6) the gluon mass acts as an 
IR regulator, making the integral convergent. Adding the three contributions together one 
obtains the black continuous curve of Fig. 5, yielding the IR behavior (3.10) with c^i = 0.33; 
this value should be compared to the value Cgog = 0.15 obtained from the SDE analysis. 
Given that the two values are derived from two a priori completely distinct methods, we 
End the proximity between the two values rather encouraging. 
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FIG. 6: One-loop diagrams contributing to the auxiliary function T appearing in the gluon 
two-point Nielsen identity (3.12). 

C. Gluon propagator 

The NI formalism may be extended in a straightforward way to the case of the gluon 
propagator. Specihcally, the corresponding NI for the gluon two-point function T^^ can be 
derived by differentiating Eq. (A20) with respect to two gluon helds, and setting afterwards 
all helds to zero. In particular, one obtains the equation 

d^^A-Ab{q) = -i'rA-xA^iq^^^A-Ab{q) -A^xAriq^^^-q)'^A-A-{q)■ ( 3 . 12 ) 

Given that T^^ is transverse to all orders, with its tree-level value given 

by r® .i,(g) = Pny{q) (see Appendix A), this identity can be further simplihed to 

read 

d^TAA{,q^) = -2irAxyi‘(g,0,-g)rAA(g^), (3.13) 

where the color structure has been factored out, and we have dehned 

rAxA*(g,0,-g) = ^^P'^''(g)rA^xAj(g,0,-g). (3.14) 

One can appreciate how the above identity works by evaluating it at lowest order in 
perturbation theory. The diagrams contributing to the function TaxA* at the one-loop level 
are shown in Fig. 6; then the rhs of Eq. (3.13) reads 

0> -iWaIaAi) = g^CAq^P'iq) |^ ~ / kA{k+ qf 

To complete the comparison, note that n^i^(( 3 ') has been evaluated in [46] [see Eq. (2.56)]; 
its derivative with respect to ^ coincides with the result (3.15), once we take into account 
that = -TA^AAq)- 
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Next, we consider the ^ 1 limit of Eq. (3.13), obtaining 

=-2ir(i^^*(g,0,-g)A“^(g2); T AA{q^) = (3.16) 

The rhs of this equation can then be evaluated within the one-loop dressed approximation, 
yielding the expression 

One notices again the presence of a massless ghost loop, which implies in turn the divergent 
IR behavior 

'9cA"^(g^)|, 0 ~ '9^m2(g^)| ~ c^ilog^ x 771^(0), (3.18) 

where the first expression on the rhs originates from the fact that when -A- 0, 
A“^(g^) —)■ w?{q^) [see Eq. (2.5)], and ml{q^) denotes the dynamical gluon mass in the 
Landau gauge [11-13]. 

The appearance of this particular behavior may be indeed conhrmed numerically. Specif¬ 
ically, after passing to the Euclidean metric and introducing spherical coordinates, one 
obtains 

„2(~i PIT POO 1 

-2iry,.(5,0,-5) I 

+ TA / d» (3 - sin^») siri^» / dy ,, ^T ,.(y)FL(zi) 
oyz^T^j Jo Jo R 

= {d) + (e), (3.19) 

which can be evaluated using the Landau gauge propagator and ghost dressing function 
introduced before. The results are shown in Fig. 7; one observes a logarithmic IR divergence 
in diagram (d), while the IR finiteness of diagram (e) is due to the presence of the dynamical 
gluon mass. When summing everything together the IR behavior is indeed the one described 
by Eq. (3.10), with Cni = 0.13. 

Finally, it is rather interesting to consider how the IR divergence found in (3.18) might be 
reconciled with the underlying assumption of an IR hnite gluon propagator. Given that, at 
present, the dynamical equation that describes the gluon mass has only been derived in the 
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FIG. 7; (color online). Contributions of the one-loop dressed auxiliary functions to the gluon 
two-point function Nielsen identity. Also in the gluon case the IR region is perfectly described by 
the predicted 61ogg^//U^ behavior, now with Cni = 0.13. 


Landau gauge®, one may only proceed by postulating an Ansatz for rn?{q^) that would satisfy 
(3.18), and study its consequences at the level of the corresponding gluon propagators. 

One such possibility is given by the following Ansatz for the .^-dependent mass function® 

,2 X'? 


m\q^) = 


a(0 + c(0 


log 


1 . 

/i2 


m^{q 


(3.20) 


with 


— Oo + + • • • ; — CiC + • • • • (3-21) 


Notice that the (resummed) behavior ~ (q'^/r^)^ has been also observed when studying the 
gfp-dependence of fermion propagators through LKF transformations [40, 51]. 


® For related studies in the Coulomb gauge, see [47-50] 
® A simpler Ansatz would have been 


m^{q^) = 


«(0 + c ( 0 (^ 


.,2 \ ? 




with 

a(0 = 1 - Co - H-; c(0 = co H- 

and Co = Cni ~ 0.13. In this case, however, the limits ^ 0 and —>■ 0 do not commute, contrary to 

what happens with the Ansatz (3.20). 
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FIG. 8: (color online), ^-dependence of the gluon mass (left panels) and gluon propagator (right 
panels) as predicted by the Ansatz (3.20) for ai = —0.2 (upper panels) and ai = 0.2 (lower panels). 

Evidently, choosing oq = 1 and Ci = Cni = 0.13 ensnres that 

~ (l + aiOmL(0), 

^^~^clog^ x (3.22) 

in agreement with (3.18); in addition, small valnes of Oi would make the and Landau- 
gauge propagators and dynamical masses to be rather close to each other, justifying in 
retrospect our replacing A by Al when solving the ghost SDE. 

Evidently, within this approach, the sign of the coefficient oi remains undetermined. 
This sign, in turn, controls the leading behavior of the gfp-dependence of the gluon mass 
(and correspondingly of the propagator) in the deep IR: a positive Oi implies an increasing 
(decreasing) mass (propagator), while for oi negative the behavior would be reversed. This 
is shown in Fig. 8, where the left panels depict the ^-dependence of the gluon dynamical 
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FIG. 9: (color online), ^-dependence of the gluon mass (left) and gluon propagator (right) as 
predicted by the Ansatz (3.20) for m = 0. 

mass (3.20) for the two values ai = 0.2 (upper-left) and ai = —0.2 (lower left), while the 
corresponding gluon propagators, obtained from the relation A“^(g^) = 
are shown on the right panels of the same hgure. 

Notice that the case Oi = 0 would be particularly interesting, as it would imply that, 
at leading order in the gluon mass and propagator coincide in the IR with the cor¬ 
responding quantities computed in the Landau gauge. In addition, as can be appreciated 
from Fig. 9, the .^-dependence over the entire range of momenta would be minimal. 

IV. CONCLUSIONS 

In the present work we have analyzed the nonperturbative behavior of Yang-Mills Green’s 
functions quantized in a linear covariant gauge, paying particular attention to its dependence 
on the parameter ^ characterizing this class of gauges. We have hrst focussed on the ghost 
two-point function and shown that, within a well-dehned set of approximations, the solutions 
of the corresponding SDE for ^ > 0 are such that the dressing function F{q^) vanishes as 
q^ —)■ 0; this is in sharp contrast to the Landau gauge case (^ = 0) where F{q^) is known to 
saturate in the low momentum region. The particular IR behavior found for F{q^) turned 
out to be in notable agreement with that obtained from the Nielsen identity satished by this 
function, within the one-loop dressed approximation and for 1. The NI analysis has 
been then extended to the gluon two-point function, and shown to predict the same kind 
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of logarithmic divergence for the derivative with respect to ^ of the dynamical gluon mass, 
~ clogg^/p^ x ml{0). A particular example of a rn?{q^) that reconciles this 
behavior with the assumed saturation of the gluon propagator away from the Landau gauge 
was given, and its main features were studied numerically. 

Undoubtedly, lattice simulations would be crucial for verifying or amending the hndings 
of this preliminary SDE study. As already mentioned, exploratory simulations in the linear 
gauges have already been carried out for the gluon propagator [34]; it would be interesting to 
extend them to larger values of in order to determine whether the observed IR saturation 
persists. Furthermore, the IR suppression of the ghost dressing function predicted here may 
serve as a dehnite reference when attempting to simulate the ghost sector of the theory. 

From the point of view of the SDFs, one may envisage various improvements. To begin 
with, the replacement of the fully-dressed ghost-gluon vertex by the tree-level expression 
inside the ghost SDF ought to be ameliorated. This, in turn, would require the treatment 
of the corresponding vertex SDF, for a general in the spirit of the analysis presented in 
the Landau gauge [43]. To be sure, subtractive instead of multiplicative renormalizability 
is another longstanding drawback in practically all types of SDF analysis; however, given 
that this problem cannot be even solved within the context of the (easier) Landau gauge, 
the prospects for a notable rehnement in this particular direction seem rather reduced. 

It is also clear that additional theoretical work at the level of the gluon propagator is 
an absolute requirement before any hrm statements could be made. In particular, no study 
related to the possibility of gluon mass generation away from the Landau gauge has been 
carried out to date; in the present work we have simply assumed the realization of this 
scenario, based almost exclusively on the limited lattice evidence of [34]. In particular, it 
would be essential to derive the dynamical equation that governs the evolution of the gluon 
mass for an arbitrary and explore the type of solutions it might admit. This task is 
technically rather complex, mainly due to the proliferation of terms with respect to the 
Landau gauge case. Calculations in this direction are already in progress, and we hope to 
report progress in the near future. 
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Appendix A: Derivation of the Nielsen identities 


The action of the SU(A) Yang-Mills theory can be written as the sum of three terms, 


p(o) 


*S'ym + St 


GF+FPG 




BV? 


(Al) 


where the hrst term corresponds to the classical action 

Sy„ = - J /d** KyFr-, F;, = d,A: - a„Al + gr’’°A‘'^Al, (A2) 

while the second to the gauge hxing and its associated Faddeev-Popov action, written as 

Agp+ppg = sjd^x c“ - ^6“^ . (A3) 

In the equation above is the Nakanishy-Lautrup multiplier, c“ the antighost held, while 
jra represents, for the moment, an arbitrary gauge hxing function. The only restriction on 
this latter function is that it allows for the inversion of the tree-level two-point functions of 
the A-b sector, thus yielding the held propagators (in what follows we will use an oh-shell 
formalism, keeping explicitly the b helds, which, otherwise, can be eliminated by making use 
of their trivial equation of motion). Finally, s is the BRST symmetry operator that acts on 
the elementary helds according to 

sc“ = &“; = 0, (A4) 

'- v -' 2 

pab 

with the usual covariant derivative. 

As can be explicitly seen above, the BRST variations of the gauge and ghost helds are 
non-linear in the quantum helds; their renormalization is ensured by the introduction of 
external sources, known as antihelds, in the third term of (Al), reading 

Sbv = j dV + c>c“). (A5) 
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The tree-level action (Al) will then satisfy the Slavnov-Taylor (ST) identity 


5(r®) = 0; 


<5r(0)5rw 5rW5r(0) 

H-h & ' 


6A^ 


5c* (5c° 


(5c“ / 


(A6) 


As the theory is anomaly-free, the ST identity (A6) holds also for the full vertex functional T. 

If we extend the BRST to include also the gauge parameter we obtain an extended ST 
identity that gives control over the gfp-dependence of the Green’s function of the theory [36, 
37]. Writing’^ 

= x; sx = 0, (A7) 


one obtains that 

^GP+PPG = I d"a; - hi - + j d^x c“ , (A8) 

and therefore the tree-level action satishes the extended ST identity 

5',r.».)^o; 5',r.».)^5(r.».)Ax^. ,A9) 


Again, the identity above is valid for the full vertex functional T; taking then a derivative 
with respect to y and setting it to zero afterwards, one obtains the NI 


dV 


/dV ( 

6^T 

6T 

6^T dr 

dr d^T 

6“ I 



x=0 

/ \5A*-dx5A^ 

dxSA^^ 

5A^ 

dx5cl 5c°- 

dc* (9ydc“ 

dx5(A J 

x=o 


(AlO) 


1. Linear covariant ganges 

Even though the NI (AlO) holds irrespectively of the gauge fixing functional chosen, we 
will specialize from now on to the case of linear covariant (or R^) gauges, which are identihed 
by the choice of the following gauge hxing function 

jra ^ 

and, in our conventions, the non negativity condition on ^ [32], needed to ensure that the 
(Euclidean) path integral over the b helds is Gaussian. 

^ A pair of variables (m, v) such that su = v and sr = 0 is called a BRST doublet; notice that also c and b 
form such a doublet. 
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iTgi,ba^{-q,q,0) = 



FIG. 10; Feynman rules for the 6-sector. Notice that, due to the 6 equation (A14), there are no 
possible quantum correction to these Feynman rules. 


Thus in the two-point gluon sector the gauge fixing yields the tree-level propagators 
AA, Ab, bb given respectively by 

= iA“'’ = 0. (A12) 

p L ^ 

For the ghost sector the tree-level propagator is instead written as 

iD‘^\q) = (A13) 

q2 

Now observe that the 6-equation 

A = a^Al - {6“ + (A14) 

implies that the 6 dependence is confined at tree-level, and so will the mixed bA propaga¬ 
tor (and any vertex involving the 6 field for that matter). Thus, beyond tree-level the 

only non-trivial propagators will be 

iAJ(«) = \p,^{q)A{q‘) + ^; tD‘\q) = AA, (A15) 

L g/ J qZ 

where F{q^) is the ghost dressing function. The Feynman rules for vertices involving fields 
and/or antifields can be found in [44]; they need to be supplemented with one more rule, 
describing the coupling of the y source to a 6 and a c fields, which can be read off directly 
from Eq. (A14): 

iFg6f,a^(-g, g, 0) = (A16) 

As already mentioned, this vertex will not receive quantum corrections, and will be com¬ 
pletely fixed by its tree-level value given above. The Feynman rules involving the Nakanishy- 
Lautrup multiplier 6 are summarized in Fig. 10. 

In addition, the Faddeev-Popov equation 


6F 6F 1 

(5c“ 5A*^ ~ 2^ ’ 


(A17) 
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implies that beyond tree-level the dependence on c of the vertex functional T can only be 
realized through the combination A*^ = A*^ + indeed, if F = r[yl*] one has 


ST 


<5r SA^(y) 


6r 


5(A{x] 


= /Attt— dl5{x-y) = -d- 


5Al{y) 5c-{x) 


^Al{y) 


6r 




^6A^^{x) 


(A18) 


with the last step due to the linearity of the held transformation employed. Making then 
the change of variable A* —)■ A*, and introducing the reduced functional F through 


F 



b^d^^A 


a 



(A19) 


one can restrict the sum over helds appearing in the rhs of Eq. (AlO) to the pairs (A“, A*“) 
and (c“, c*) alone. In the NI analysis carried out in this work we have use only ‘tilded’ 
quantities, and therefore suppressed this symbol everywhere. Incidentally, notice that it 
is Eq. (A19) that implies the tree-level result F® .(,(g) = iq^5‘^^Pnu{q). 

The hnal form of the NI used is then written as 


dV 


= / d'^a; 


5F (J^F (J^F (5F (J^F 5F (5F (5F 
dxbAa dxbA*^ 5Aa dxScl be* dxbc° 


(A20) 


Using the technique developed in [52], one can write the complete solution to the NI 
above [53]. Rewriting Eq. (A20) as® 


dV 


= / (Tx 


x=o 


(5T (5F (5T (5F (5T (5F (5T (5F 

I r .. r ... r ... r 


5A*^5A'i, SAfbAT Sc^^Sc: Sc* Sc 

11 fi 


X=0 


its full solution is given by [53] 


dT 

(A21) 


r = VU”r„; r„=[Ajr 

^ nl 


oJI^=o: 


n>0 


(A22) 


where Fq = F|^^q is the vertex functional in the Landau gauge, and in the gauges the 
Lie operator A,i, reads 


A>i,X 



/ SX (5T SX (5T (5A(5T (5A(5T\ dX 


(A23) 


® The sign differences with respect to [53] are due to the different conventions used. In particular, our 
Yang-Mills action is obtained from the one of [53] through the replacements: c —>■ —c, b —>■ —b, c* —>■ —c* 
and a —>■ —^ (which also implies 9 —>■ —y when introducing the doublet partner of the gfp parameter). 
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If 1 one can linearize Eq. (A22); the coefficient of the linear term Ei is then obtained 
by applying the Lie operator on the Landau vertex functional Lq. As the latter does not 
depend on the gfp within the linear approximation one has 


= El, 


(A24) 


with 


El = / d^x 


= / d^x 


(5Eo 5^ (5Eo 5^ (5Eo 5^ (5Eo 5^ 

I r ^_ r \ U, I r . 1 . r ^ I 


SA- SAT 


STr 


(J^Er 


SA*- SAH 
S‘^U 


Scl Sc°- 


5c°- Sc* 


5En 


^^En STr 


«=0 

^ho <5Eo 


(A25) 


SA*^ dxSAa dxSA*^ SAa dxScl Sc^ Scl dxScf 
We then see that the approximation employed in this work on the full NIs are equivalent to 
differentiating Eq. (A24) with respect to a ghost and an antighost [Eq. (3.8)], or two gluon 
helds [Eq. (3.16)]. 
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